Split vortices in optically coupled Bose-Einstein condensates 
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We study a rotating two-component Bose-Einstein condensate in which an optically induced 
Josephson coupling allows for population transfer between the two species. In a regime where 
separation of species is favored, the ground state of the rotating system displays domain walls with 
velocity fields normal to them. Such a configuration looks like a vortex split into two halves, with 
atoms circulating around the vortex and changing their internal state as they cross the domain wall. 
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Vortex formation is generally viewed as an unequivocal 
signature of superfluid motion in atomic Bose-Einstein 
condensates. Both in one- and two-component systems, 
a number of vortex-like structures have been created 
and observed ^ that confirm theoretical predictions 
1^, ^, ||] . The Josephson effect between two weakly cou- 
pled condensates is another paradigm of superfluid trans- 
port that so far has received less attention [|[ !?[ ^. In 
this article, we present a combined study of these two 
fundamental properties of superfluid systems: vortex and 
Josephson dynamics. We study the ground state proper- 
ties of a rotating double condensate system in which the 
combined role of vortex formation and optical coupling 
between two different hyperfine states gives rise to a rich 
physical behavior. 

A crucial consequence of the internal Josephson cou- 
pling is the generation of an effective attraction between 
both atomic species due to the energy that atoms gain 
by choosing a symmetric mixture of the coupled internal 
states. Therefore, the most interesting physics is reached 
by combining this coupling with setups that otherwise 
favor species separation, such as a particular combina- 
tion of atomic scattering lengths H, [l0| o r a separation 
of the respective confining potentials In this type 

of setups we find that the effective attraction due to the 
optical coupling causes an increase in the thickness of 
the domain wall where the two components physically 
overlap. More important is the fact that, if we add ro- 
tation to a Josephson coupled condensate with separate 
domains, atoms can use the domain wall to mutate their 
internal state and shift between components in a conti- 
nous way. Combining this persistent current in the inner 
space with a persistent current in real space, the double 
condensate may now easily create a vortex core for each 
component in the region where that component has a low 
density. Even for otherwise small values of the Joseph- 
son coupling, these split vortices support a net mass flow 
comparable to that of conventional vortices. The forma- 



tion of these novel structures is less costly in terms of 
angular speed because the vortex of a given component 
is formed not within its own domain but in the opposite 
one, where its superfluid density is low and the cost in 
kinetic energy is therefore small. 

To study this type of vortex structure, we analyze first 
a rotating two-component condensate without optical 
coupling (i.e. with impenetrable domain walls) and show 
that its behavior is essentially that of a one-component 
system with a displaced axis of rotation. Then we show 
that vortex formation is strongly inhibited because of 
the ability of the system to gain angular momentum by 
merely distancing itself from the rotation axis. The pic- 
ture changes qualitatively when a Raman coupling is in- 
troduced to permit coherent hopping between the two 
internal states. Because the flow of particles in a given 
state is no longer a conserved quantity, the circulation 
lines of a component can cross the domain wall with a 
concomitant decrease in their supporting superfluid den- 
sity. This results in a global structure of two asymmetri- 
cal vortices where matter is efficiently transported both 
in real space and within the hyperfine doublet. 

The model.- In this paper we focus on double conden- 
sate systems such as those made of Rb in JILA [Q , but 
this time in rotating traps and with a permanent optical 
coupling between the species. In the rotating frame of 
reference Q the Gross-Pitaevskii equations for the con- 
densate wavefunctions read 



ihdr'^i = 



*i - i/if7/?*2, (la) 



H2+Y.,U23\^j? *2 - \hQR-^i. (lb) 



Here Hj correspond to the non-interacting Hamiltoni- 



ans Hi ■ 
V{v) = 



2m 



V{r — ri_2) — hflLz + Mi^2, where 
is the trapping potential, is the cen- 
ter of the trap for component j, and fl is the angular 
speed. The mutual interactions Uij = ATrh?aij /m, are 
proportional to the s-wave scattering lengths a^ . The 
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FIG. 1: (a) Density profiles in scenario A for f2 = 0. We show 
a grayscale plot of IV'il'^i (dark), and a contour plot of \1p2f, 
(solid lines). The trap separation is indicated with dashed 
lines, (b) Type A condensates for = 0.135. Grayscale 
plot for the phase {9 — axgipi), arrows for the direction of 
the velocity field (v — VS), and contour lines for ■ (c) 
Separation of the atomic clouds and (d) angular momentum 
as a function of the angular speed for scenarios A (circles, 
solid line) and B (crosses, dashed line). All magnitudes have 
been adimensionalized. 



last term in each equation models an optical coupling 
between species. This is a Josephson-type coupling that 
allows atoms to change their internal state. Since we 
focus on stationary configurations, the normalization of 
each wavefunction is fixed, Ni — J \'i>i{r)\^dr. Finally, 
hSi is a tunable energy splitting in the hyperfine space, 
which controls the population of each component. 

In this paper we restrict our attention to axially sym- 
metric two-dimensional configurations which may de- 
scribe pancake type traps p^ . In such traps the non- 
linear parameter is affected by a corrective factor |l^ 
and A'^i and N2 are interpreted as effective number of 
particles. To ease the analysis we introduce dimension- 
less variables x = r/oo and t = t/T . With this change 
and '4>j{'X-, t) — ' ^'j(r, r), we write the equations for 
stationary states 



-01 - AV'2, 



(2a) 
(2b) 



with rescaled Hamiltonians 



ni,2 = i [-A + (x - xu2f + y^] + in {xdy ~ y9^) . 



(3) 



FIG. 2: Scheme of an ofT-axis rotating trap. Both the center 
of the trap and the trap itself rotate with angular speed Q. 



In Eq. (H) the splittings 5i are included in the chemi- 
cal potentials. The parameter A = 217^/^ measures the 
intensity of the optical coupling. Although it is widely 
tunable, in our work it will be spatially uniform and at 
most of order unity. We have solved numerically Eq. (H), 
looking for the solutions that have lower energy, the so 
called ground states. Each of such solutions represent a 
stable, experimentally realizable configuration. 

We will consider two different scenarios in which the 
double condensate exhibits domain walls, and which give 
qualitatively similar results. The first case, which we call 
"setup B", corresponds to a situation in which gu — 
<?22 = N,gi2 = 2N, with a choice of iV = 38. In this case 
the inequality > 5ii522 is satisfied and the domains 
form spontaneously pO| with no need for trap separation 
nor splitting, i. e. xi,2 = 0, fii — fi2- 

The second and most important scenario, which we 
call "setup A" corresponds to the case of ^^Rb |l[] with 
parameter values ( g^l gll) = { 0.94 o!97 ) ^ (i-^- ^1 — 
N2 = N), and a typical effective value of = 100. This 
type of condensates has been studied in many previous 
experimental and theoretical works. For ^^Rb, the in- 
equality mentioned above is close to saturation and 
a separation xi = —X2 — 1 ensures the formation of two 
different domains. The external splitting takes small val- 
ues, fJ-i ~ fJ-2 ^ 0.01, and it does not influence the results. 

Rotation without Josephson coupling.- Here we con- 
sider the case A = 0. The existence of domain walls im- 
plies that the motion of species is spatially constrained. 
Therefore, if we impose some angular speed to the traps 
containing the condensates, each cloud tends to slip tan- 
gentially to the domain wall. The density distribution is 
similar to the case without rotation, but now, due to the 
centrifugal force, the species separate a little and gain lin- 
ear speed. Together with the deformation of the clouds, 
this mechanism permits the acquisition of a large amount 
of angular momentum without generating vortices [Fig. 
|l|(d)], which only form at very high fi. 

In Fig. ^a) wc show the ground state of the rotating 
double condensate system for A = 0, showing a structure 
of domain walls that persists in the presence of rotation 
and prevents the formation of vortices, as revealed by the 
phase and circulation patterns of Fig. |l](b). From Fig. 
p](c) it is evident that the both the mutual repulsion and 
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the rotation of the trap increase the separation among 
species. The flow pattern presents a curvature which is 
too small [Fig. |l|(b)], and the gain of angular momen- 
tum [Fig. |l|(d)] is due instead to the displacement of the 
condensate away from the origin. 

To get a deeper understanding of this configuration it is 
useful to study the off-axis rotation of a single component 
condensate. It experiences a potential V{x,t) = i(x — 
ro{t))A{t){x—ro{t)) with A{t) and ro{t) rotating at speed 
ri around a displaced axis, as depicted in Fig. ^j. Thanks 
to a particular symmetry of the nonlinear Schrodinger 
equation with harmonic potential |l4| , any solution can 
be written as 



V'(x, t) = 0(x - u, i)e™"-V'^+/(*) 



(4) 



where (j)(x — u) is centered on the center of mass u and 
satisfies the Gross-Pitaevskii equation (j^ for ro = 0. In 
the case of radially symmetric traps, A = ui^, we find 
(1 — fl^/tt'^)u = Tq. This tells us that the behavior of 
a condensate in an off-axis rotating trap is qualitatively 
similar to that of a condensate in a centered trap: Vor- 
tices nucleate at similar angular speeds and they all ap- 
pear inside the condensed cloud. The difference is that 
now the condensate has an additional source of angular 
momentum due to its displacement with respect to the 
origin, cc fl\u\'^. 

These arguments, which are rigorous for the single con- 
densate system, may be extended to the two-component 
condensate case where the overlap between clouds is 
small, meaning that vortices should appear in the center 
of the corresponding clouds even when the trap centers 
are displaced with respect to the axis of rotation. 

Role of Josephson coupling.- In order to explain the 
role of Josephson coupling let us first consider situation 
A. It is easy to show that the Josephson terms favor 
energetically the mixing of both species. In the limit of 
strong coupling, they coexist in space even when their 
traps are separated. This is most intuitively appreciated 
by inspecting the energy functional associated to Eq. (0) 



£^[^1,^2] 



[1pl^-^l^pl + ip2'H2'4'2\ 



ARe (ipi'4'2 



(5) 



The coupling term by itself is minimized with a solution 
such that argipi = arg-02. The analogy with the off-axis 
rotation of a single condensate [see Eq. (^] suggests the 
variational wavefunction 

M^,y) « e-(^+«)V2-yV2g-»^a. (6b) 
Substituting this ansatz into Eq. (S) one gets an effective 
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FIG. 3: Rotating two-component condensate with coupling 
A = 0.1, scenario B with TV = 100 and angular speed Q — 0.2. 
(a) Contour lines for the moduli of the wavefunctions \'4^i\ 
(solid line) and |'i/'2| (dashed line), (b) Contour levels for 
\tpi\, phase {9 = argf/ji, grayscale), and velocity field (vi = 
Vargi/)i, triangles), (c) Condensate fluxes jl*' = 
together with intensity of the optical coupling, A|^/)i 11/12 j, along 
the line y = —0.097. (d) Position of the vortex for a type B 
condensate with = 0.28, as a function of the coupling. All 
units have been adimensionalized. 



energy 



£: - i(u - d)2 + + nuv - Ae"" + .giaCe 



(7) 



where C is of order unity. Minimization with respect 
to the velocity v leads to v — —^lu, and thus fl favors 
separation. More importantly, it is clear from (^ that A 
effectively decreases the repulsion between condensates. 
Thus, the separation u decreases with the strength of the 
optical coupling and becomes zero for a strong enough 
value of the Josephson couphng, Ac. Both the repulsive 
interaction among bosons and the rotation of the trap, 
tend to inhibit mixture of species, thus increasing the 
value of Ac. 

Vortex formation.- The structure of matter flow in 
the condensate suffers a drastic trasformation when the 
Joscpson coupling is introduced. From numerical solu- 
tions of Eqs. ^ for setups A and B, we see that any 
nonzero coupling allows the formation of vortices [Figs, 
^(a) and^. These vortices involve a matter flow which is 
orthogonal to the wall separating both condensates [Fig. 
^(b)]. We note that the exact location of both vortex 
cores depends on the intensity of the coupling [Fig. ||(d)]. 

The mathematical reason for this striking change is 
that, with a finite value of A, the matter flow of each 
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angular speed. The cliffs on the surface are due to the 
nucleation of successive vortices. The separation of these 
vortices from the a; = domain wall decreases very fast 
with increasing A: Fig. ||(d) shows that, already for A = 
0.05, the vortices become visible. As A — > 0+ for fixed 
f2, vortices disappear continuously moving their core to 
infinity. To create vortices at A = 0, one needs to reach 
a much higher angular speed, fie, comparable to that of 
single component condensates Q . The upshot is that the 
optical coupling permits the formation of split vortices 
with important mass flow at lower angular velocities. 

Internal Josephson dynamics.- The Hamiltonian (5) 
may be written as that of a nonrigid pendulum ||] with 
an effective interaction energy Ec = X^ij 5u '''^ ("1) > 
where ct^ = / ji/'iPlV'j and an effective Rabi frequency 
lor = njisi2, with si2 = /'0iV'2- The conclusion is that, 
for the setups considered here (A ^ 0.5), the internal 
two-state dynamics lies in the collective Josephson regime 
(2ujr/N <^ Ec ^ Nujr/2), while in the JILA experiment 
Ig , where no walls are formed, the same dynamics lies in 
the noninteracting Rabi limit {Ec <C 2ujj^/N). 

This work has been partially supported by Ministerio 
de Ciencia y Tecnologia under grants BFM2000-0521 and 
PB96-0080-C02. 



FIG. 4: (a-b) Angular momentum per particle, h, versus Q 
and A. (a) type A condensates (b) type B condensates, (c-e) 
Density plots and (f-h) phase plots of states of a type A rotat- 
ing condensate: A = 0.315, O = 0.316 (c,f), = 0.447 (d,g) 
and Q = 0.474 (e,h). In (c-h) we plot both species separately 
(although they overlap). AH magnitudes are dimensionless. 



species is no longer conserved. Assuming a divergence 
free flow, the equation of continuity along the trajectory 
of a boson becomes 



dl 



-vi = -ARe(V'i'02) = - 



dl 



-V2. 



(8) 



Thus, as the current line is closed around the origin, there 
is an exchange of bosons among components. A typical 
boson flowing around the origin suffers a transformation 
from state |1) to |2) and viceversa as it completes a circle. 

Physically, there is a fundamental reason why bosons 
are transferred from one component to the other. For 
A = 0, vortices must lay inside each condensed cloud, 
due to current conservation. However, for any nonzero 
coupling, vortices may appear outside the bulk of the 
clouds, at a variable distance from the origin. Placed 
at low density regions, the twist of the phase requires 
less kinetic energy J argipi)"^ and provides more 

angular momentum than the mere separation of clouds. 

In Figure ^a-b) we plot the angular momentum per 
particle as a function of the Josephson coupling and the 



[1] 
[2] 

[3] 
[4] 

[5] 
[6] 
[7] 
[8] 



[9] 
[10] 

[11] 

[12] 
[13] 



M. R. Matthews, B. P. Anderson, P. C. Haljan, D. S. 
HaU, C. E. Wieman, E. A. Cornell, Phys. Rev. Lett. 83, 
2498 (1999). 

K. W. Madison, F. Chevy, W. Wohlleben, J. Dalibard, 
Phys. Rev. Lett. 84, 806 (2000); F. Chevy, K. W. Madi- 
son, J. Dalibard, Phys. Rev. Lett. 85, 2223 (2000); J. R. 
Abo-Shaeer, C. Raman, J. M. Vogels, and W. Ketterle, 
Science 292, 476 (2001). 

See e.g. A. Svidzinsky, A. L. Fetter, Jour. Phys. B 13, 
R135 (2001), and references therein. 
D. L. Feder et al. Phys. Rev. Lett. 86, 564 (2001); F. 
Dalfovo, S. Stringari, Phys. Rev. A 63, 011601 (2001). 
J. J. Garci'a-Ripoll and V. M. Perez-Garcia, Phys. Rev. 
A 64, 013602 (2001). 

D. S. Hall, M. R. Matthews, C. E. Wieman, E. A. CorneU, 
Phys. Rev. Lett. 81, 1543 (1998). 

B. P. Anderson and M. A. Kasevich, Science 282, 1686 
(1998). 

For a review, see F. Sols, in Proceedings of the Interna- 
tional School of Physics "Enrico Fermi", M. Inguscio, S. 
Stringari, and C. E. Wieman, eds. (lOS Press, Amster- 
dam, 1999), p. 453. 

B. D. Esry, C. H. Greene, Nature 392, 434 (1998). 

S. Coen, M. Haeherman, Phys. Rev. Lett. 87, 140401 

(2001). 

D. S. HaU, M. R. Matthews, J. R. Ensher, C. E. Wieman, 

E. A. CorneU, P hys. Rev. Lett. 81, 1539 ( 1998). 



arXiv: condinat/0104549 



A. Gorlitz et al., 

J. J. Garci'a-RipoU, G. Molina-Terriza, V. M. Perez- 
Garci'a, and L. Torner, Phys. Rev. Lett. 87, 140403 
(2001). 

[14] J. J. Garci'a-Ripoll, V. M. Perez-Garcfa, and V. Veksler- 



chik, Phys. Re v. E. in press 
M. Linn et al. 



arXiv : condmat/0103489 



